function [points, weights] = gauss_lengdre_loc(V, int_type)
Jacobi = [V(2,:)-V(1,:); V(3,:)-V(1,:)]';
J_det = det(Jacobi);
[gauss_points_ref, gauss_weights_ref] = gauss_lengdre_ref(int_type);
weights = abs(J_det) * gauss_weights_ref;
points(:,1) = V(1,1) + Jacobi(1,1)*gauss_points_ref(:,1) + Jacobi(1,2)*gauss_points_ref(:,2);
points(:,2) = V(1,2) + Jacobi(2,1)*gauss_points_ref(:,1) + Jacobi(2,2)*gauss_points_ref(:,2);
end